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ABSTRACT 



Aims. Determination and analysis of electron densities and filling factors of the difl[^use ionized gas (DIG) in the solar neighbourhood. 
Methods. We have combined dispersion measures and emission measures towards 38 pulsars at distances known to better than 50%, 
from which we derived the mean density in clouds, A'c, and their volume filling factor, Fy, averaged along the line of sight. The 
emission measures were corrected for absorption by dust and contributions from beyond the pulsar distance. 

Results. The scale height of the electron layer for our sample is 0.93±0. 13 kpc and the midplane electron density is 0.023 ±0.004 cm"^, 
in agreement with earlier results. The average density along the line of sight is (n^) = 0.018 ± 0.002 cm"'' and nearly constant. Since 
(We) = F^Nc, an inverse relationship between Fy and N^ is expected. We find FyiN^) = (0.01 1 ± 0.003) ;\f-i 20±o.i3^ which holds for the 
ranges N^ = 0.05 - 1 cm"^ and Fy = 0.4 - 0.01 . Near the Galactic plane the dependence of Fy on N^ is significantly stronger than away 
from the plane. Fy does not systematically change along or perpendicular to the Galactic plane, but the spread about the mean value 
of 0.08 ± 0.02 is considerable. The total pathlength through the ionized regions increases linearly to about 80 pc towards |z| = 1 kpc. 
Conclusions. Our study of Fy and A'c of the DIG is the first one based on a sample of pulsars with known distances. We confirm the 
existence of a tight, nearly inverse correlation between Fy and A'c in the DIG. The exact form of this relation depends on the regions 
in the Galaxy probed by the pulsar sample. The inverse Fy-N^ relation is consistent with a hierarchical, fractal density distribution 
in the DIG caused by turbulence. The observed near constancy of (n^) then is a signature of fractal structure in the ionized medium, 
which is most pronounced outside the thin disk. 

Key words. Galaxy: disk - HII - ISM: clouds - ISM: structure 



1. Introduction 

The full-sky Ha maps of the Milky Way that have recently be- 
come available (Dickinson et al. 120031 Finkbeiner l2003 1 Bennett 
et al. l2003i: Hinshaw et al. 120071) show a spectacular variation in 
intensity and structure. The classical H II regions in the disc are 
suiTounded by extended, diffuse ionized gas (DIG) visible up 
to high latitudes, and supernova explosions, shocks and turbu- 
lence in the interstellar medium (ISM) have produced a wealth 
of cloud shapes, filaments, shells and voids. 

Madsen et al. (2006) found the variations in structure to be 
accompanied by variations in electron temperature. They derived 
a mean temprature of the DIG of about 8000 K with strong dif- 
ferences between lines of sight. 

A number of authors derived a scale height of the DIG in 
the solar neighbourhood of about 1 kpc from the increase of the 
dispersion measures of pulsars with increasing distance to the 
Galactic plane (Reyn oIds [l991bi Bhatta charia & Verbunt'1991', 
Nordgren et al.[19^ Gomez et al. 120011 Cordes & Lazio 2003). 
Using dispersion measures towards 4 pulsars in globular clus- 
ters, Reynolds ( |1991al ) obtained a filling fraction of the DIG of 
S 0.2 through the full layer and a mean density in clouds of 
about 0.08 cm"^. This was recently confirmed by Berkhuijsen 
et al. ( 120061 ) from a much larger pulsar sample. These authors 
also found an increase of the filling factor with distance from the 
plane. Already the scanty data that were available to Kulkarni & 
Heiles ( 119881 ) had indicated such an increase. 

Turbulence causing hierarchical, fractal structure leads to an 
inverse relationship between the mean filling factor along the 
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line of sight, Fy, and the mean density in ionized regions, A^c 
(Fleck [T9961 Elmegreen [19981 [19991 ). Pynzar ( [T9931 ), who was 
the first to investigate this relationship, derived Fy oc N^'^''' for 
a combination of DIG and H II regions, which was confirmed 
by Berkhuijsen ( 119981 ). In an extensive study, Berkhuijsen et al. 
([20061 hereafter called BMM) obtained Fy oc a?-i 07±o.03 f^j. jjjg 
DIG from a sample of 157 pulsars at latitudes \b\> 5°. They used 
dispersion measures from the catalogue of Hobbs & Manchester 
(2003 ; see also Manchester et al. 2005 ), emission measures from 
the WHAM survey (Haffner et al. 120031 ), corrected for absorp- 
tion (Diplas & Savage 19941 as well as for contributions from 
beyond the pulsars, and pulsar distances from the model of the 
electron distribution in the Galaxy of Cordes & Lazio ( 20021 ). 
The statistical error in pulsar distances < 3 kpc, where 75% of 
the pulsars in the sample are located, is about 25% which is 
much smaller than the intrinsic spread in the dispersion mea- 
sures and emission measures. Therefore, random eiTors in the 
model distances will not have influenced the statistical results. 
However, the distances from the model may also contain a sys- 
tematic error. Based on these distances, the radial distribution of 
pulsars has a maximum at /? = 3.5 kpc from the Galactic centre 
(Lorimer 2004; Yusifov & Kiiciik [20041 ), whereas other popula- 
tion I objects (H2, HI, HII regions) peak near R - 5 kpc. The 
smoothness of the electron density model may lead to pulsar dis- 
tances that are too large. 

So it is important to check the results of BMM for pulsars 
with distances derived from observations. Since BMM did their 
work, the sample of pulsars with measured distances has in- 
creased considerably. Our new analysis is based on 38 pulsars, 
the distances of which are known to better than 50%. A further 



Elly M. Berkhuijsen and Peter Miiller: Densities and filling factors of the DIG in the Solar neighbourhood 



decrease of the constraint on the distance errror would make the 
sample too small for a statistical analysis. We also used an im- 
proved absorption correction to the emission measure of each 
pulsar instead of the statistical correction applied by BMM. 

2. Basics and data 

2.1. Basic relations 

The expressions for dispersion measure, DM, and emission mea- 
sure, EM, towards a pulsar at distance D (in pc) can be written 
in various ways: 



DM 



EM 



Jo 
Jo 



{l)dl = <«e>£> = N^FyD = A^c^e 



{l)dl ■ 



{nl)D 



nIf, 



D 



NlU 



(1) 



(2) 



where n^il) (in cm"^) is the local electron density at point I 
along the line of sight, Le (in pc) the total pathlength through 
the regions containing free electrons (clouds, clumps) and N^ (in 
cm"^) the average density in these regions; between clouds the 
electron density is assumed to be negligible (see Fig. 1 in BMM). 
Furthermore, (ng) and (n^) are averages along D, and Fy - L^/D 
is the fraction of the line of sight in clouds, which approximates 
the volume filling factor F^ if there are several clouds along the 
line of sight (BMM). Note that the third equality in eq. (|2]i is 
only valid when the average density of every cloud n^ along the 
line of sight is the same: then (n^) = {nl}Fy - N^Fy. Thus Nc 
and Fy are approximations of the true average density in clouds 
and their filling factor 

Combining eqs. ([TI3 and (|2]l we have 



N^ = EM/DM , 
Fv = DMVEM£). 



(3) 
(4) 



The derived quantities Nc and Fy are connected by the simple 
relation (ng) = F^Nc following from the third equality in eq. ([T]). 
As the mean electron density in clouds will not be constant, our 
estimate of A^c will be too high because of the nl dependence 
of EM. Recently, Hill et al. ( 120081 1 found that the ratio between 
the most probable density and A^c varies between 0.95 and 0.4 
in the mildly supersonic cases of their MHD simulations (scaled 
to observations). Thus, our values of Nc may be too high by a 
factor ^ 2; the corresponding values of Fy will then be too low 
by the same factor. 

2.2. The data 

We collected a sample of 78 pulsars with distances known to 
better than 50% from the literature. This error is negligible com- 
pared to the large intrinsic spread in DM (factor of 2) and EM 
(factor of 4). The sample consists of 52 pulsars from the list of 
Gomez et al. (2001 ), where we used more recent distances when 
available, and 26 pulsars with parallactic distances taken from 
Hobbs et al. ( 120051 and references therein). 

The dispersion measures came from the catalogue of 
Manchester et al. (2005). To find the emission measures in the 
pulsar directions we used the full-sky Via map of Finkbeiner 
( 120031 1 corrected for extinction as described by Dickinson et al. 



( 120031 ), who argued that only one third of the dust in the line of 
sight effectively absorbs Wa emission. The extinction-corrected 
emission measure is 

EMc = EM X 0.34 A with A = exp (2.4 E{B - V)) , 

where EM is the observed emission measure and A the correc- 
tion recommended by Finkbeiner (I2003I I. Corrections of more 
than 1 magnitude are uncertain and cannot be used. Therefore 
we had to remove 38 pulsars from the sample, most of which are 
at latitudes \b\ < 5°. We also removed 2 pulsars with high EM 
and DM (J0613-0200 and J1807-H0943) indicating H II regions 
on their lines of sight. Thus for the analysis a sample of 38 pul- 
sars is available, for 23 of which a parallactic distance is known. 
The errors in the distances to 28 pulsars are smaller than 20% 
(AD/D < 0.2) and those in the distances to the other 10 pulsars 
are less than 50% (0.2 < AD/D < 0.5). As 21 pulsars are lo- 
cated at D < 2 kpc, the sample is heavily weighted to the solar 
neighbourhood; 13 pulsars are in globular clusters at distances 
> 2 kpc. The main parameters of the 38 pulsars in the sample 
are listed in Table [1] which also gives the references to the dis- 
tances used. The sightlines to these pulsars probe the DIG. 

In the direction of nearby pulsars and of pulsars near the 
Galactic plane a significant amount of the emission measure 
originates beyond the pulsar (see Fig. 5 in BMM). Reynolds 
( 119971 ) and Haffner et al. (120031 1 have shown that the vertical 
electron distribution derived from Ha observations is fairly well 
described by an exponential. Then the emission measure up to 
the pulsar, EMp, is obtained from 



EMp sin \b\(z) = EM, sin \b\ (l - exp (-|Zp|//ie)) , 



(5) 



where Zp is the distance of the pulsar to the Galactic plane and 
hf. is the scale height of nl(z). Deviations from a smooth, plane- 
parallel layer will cause scatter in EMp. The mean of EM, sin \b\ 
for our sample is 2.8 ± 0.3 cm"^ pc, which is the extinction- 
corrected emission measure through the full layer. We derive the 
scale height h^ in the next section. 

3. Scale height of electrons 

We first estimated the scale height h^ of «e(z) using Eqs. 12 
and 13 of BMM, expecting that also for our sample the av- 
erage filling factor Fy{z) would increase with |z|. This yielded 
280 pc < /le < 510pc with a best value of 390 pc, giving the 
observed maximum values of DM sin \b\ = 22 + 2 cm"-' pc (see 
Fig- [B and EMpSin|/7| = 2.8 + 0.3cm"^pc. However, for this 
entire range of /ig, Fy(z) appeared to be essentially constant for 
our sample (see Fig.|4}l), thus also the local filling factor f{z) is 
^ constant. As is easily seen, this indicates that the scale height 
of nl(z) is about half that of ng{z). 

Following BMM, we describe the dispersion measure per- 
pendicular to the Galactic plane as 



DM sin \b\(z) 



n 
Jo 



■.(z)dz 



IZpl 



f(z)njz)dz 



= foncoh{l -exp(-|Zp|//j)) , 



(6) 



Note that we write F^ and N^ where BMM used /^, and n^ 



where we assume that f(z) - /o and the local density in clouds 
riciz) is an exponential with midplane value WcO and scale height 
h; Zp is the distance of the pulsar to the midplane. 

Similarly, the emission measure perpendicular to the 
Galactic plane is 
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Table 1. Parameters of the 38 pulsars in the final sample 



PULSAR 
Jname 



LONG 

n 



LAT 

n 



DM° 

[ cm"^ pc] 



D 



AD" REF EM, EMp 



A'c 



[pc] 



[cm pc] [cm ] 



J0437-4715 


253.40 


-42.00 


2.65 


159 


5 


1 


1.8 


0.4 


0.137 


0.122 


JOS 14+7429 


140.00 


31.62 


6.12 


433 


8 


5 


1.9 


0.7 


0.116 


0.122 


J0826+2637 


197.00 


31.70 


19.45 


360 


80 


6 


4.8 


1.6 


0.081 


0.665 


J0953+0755 


228.91 


43.70 


2.96 


262 


5 


5 


3.2 


1.0 


0.351 


0.032 


J1136+1551 


241.91 


69.19 


4.86 


350 


20 


5 


1.5 


0.8 


0.155 


0.090 


J1239+2453 


252.45 


86.54 


9.24 


850 


60 


5 


1.1 


0.9 


0.100 


0.108 


J1456-6843 


313.90 


-8.50 


8.60 


450 


60 


8 


32.4 


4.3 


0.497 


0.038 


J1932+1059 


47.38 


-3.88 


3.18 


330 


10 


5 


18.5 


0.9 


0.269 


0.036 


J2018+2839 


68.10 


-3.98 


14.18 


950 


90 


5 


35.0 


4.6 


0.323 


0.046 


J2022+2854 


68.86 


-4.67 


24.64 


2300 


800 


5 


38.5 


12.6 


0.513 


0.021 


J2022+5154 


87.86 


8.38 


22.65 


1900 


250 


5 


12.6 


5.6 


0.247 


0.048 


J0922+0638 


225.42 


36.39 


27.27 


1210 


90 


7 


7.3 


5.7 


0.209 


0.108 


J1537+1155 


19.85 


48.30 


11.61 


1080 


150 




2.7 


2.2 


0.192 


0.056 


J17 13+0747 


28.75 


25.22 


15.99 


910 


80 




4.9 


2.8 


0.172 


0.102 


J1744-1134 


14.79 


9.18 


3.14 


470 


95 




21.4 


3.2 


1.005 


0.007 


J1909-3744 


359.73 


-19.60 


10.39 


1140 


50 




8.1 


4.5 


0.431 


0.021 


J0659+1414 


201.11 


8.26 


13.98 


290 


30 


4 


14.0 


1.2 


0.085 


0.567 


J2145-0750 


47.78 


-42.08 


9.00 


500 


200 


9 


1.7 


0.9 


0.098 


0.184 


J1022+1001 


231.79 


51.10 


10.25 


400 


145 




2.7 


1.3 


0.129 


0.198 


J 1024-07 19 


251.70 


40.52 


6.48 


520 


270 




2.3 


1.2 


0.183 


0.068 


J2124-3358 


10.93 


-45.44 


4.60 


250 


165 




2.8 


0.9 


0.189 


0.098 


J0030+0451 


113.14 


-57.61 


4.33 


300 


90 


2 


1.0 


0.4 


0.099 


0.146 


J1509+5531 


91.32 


52.29 


19.61 


2370 


220 


3 


1.0 


1.0 


0.049 


0.169 


J05 34+2200 


184.56 


-5.78 


56.79 


2000 


75 




23.3 


8.1 


0.143 


0.198 


J1740-5340C 


338.20 


-11.90 


71.80 


2300 


20 




37.8 


24.0 


0.335 


0.093 


J0141+6009 


129.10 


-2.10 


34.80 


2750 


150 




29.0 


5.6 


0.161 


0.079 


J1910-5959C 


336.50 


-25.60 


33.68 


4000 


250 




3.2 


3.1 


0.092 


0.091 


J 1824-245 2c 


7.80 


-5.58 


119.86 


5700 


750 




76.6 


53.1 


0.443 


0.047 


J1701-3006C 


353.60 


7.30 


114.56 


6900 


900 




59.4 


50.2 


0.438 


0.038 


J1911+0101C 


36.11 


-3.92 


202.68 


9500 


800 


12 


39.2 


29.3 


0.145 


0.147 


J1823-3021C 


2.79 


-7.91 


86.84 


8000 


800 




26.8 


24.2 


0.279 


0.039 


J1804-0735C 


20.79 


6.77 


186.32 


8400 


1900 




35.5 


31.2 


0.168 


0.132 


J1748-2021C 


7.73 


3.80 


220.40 


8400 


2100 




88.7 


61.6 


0.279 


0.094 


J1721-1936C 


4.86 


9.74 


75.70 


8600 


1100 




21.4 


20.4 


0.270 


0.033 


J0024-7204C 


305.92 


-44.89 


24.60 


4500 


250 




3.6 








J1518+0205C 


3.8 


46.80 


29.47 


7500 


400 




2.1 








J1641+3627C 


59.00 


40.91 


30.36 


7700 


380 




1.8 








J2129+1210C 


65.01 


-27.31 


67.31 


10300 


650 




5.0 









aj Hobbs & Manchester ( I2003] >: www.atnf.csiro.au/research/pulsar/psrcat 

b) Mean of absolute values of positive and negative error. 

c) In globular cluster. The last four pulsars are at |z| > 3 kpc, and only used in the determination of the 
electron scale height (see Fig.[Tll. 

References. (1) Hotan e t al. 120061 (2) Lommen et al. 120061 (3) Ch atterje e et al. |2005l (4 ) Brisk en et al. 

(5) Brisken et al. l2002l (6) Gwinn et al. ri986. (7) Chatterjee et a l.l200T] (8) Bailes et al. [T990l 

(9) Lohmer et al. l2004l (10) Stairs et al. [T998l (11) Gomez et al. 120011 (12) Heitsch & Richtler [T999l 
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■IZpl 



EJMpSinl^Kz) = r nliz)dz^ I f(z)nt(z)dz 
Jo Jo 

= /o«?o/z(l-exp(-2|Zp|//2))/2, 



(7) 



where the scale height h^ in eq. ^ is equaQ to h/2. 

We can derive the maximum of DM sin \b\ and the scale 
height h from the distribution of DM sin \b\{z) of our sample of 
38 pulsars shown in Fig.[T] A two-parameter fit yielded fon^Qh - 
2 1.7+ 1.5 cm"-' pc and /i = 0.93+0.13 kpc, giving ne(0) = /oWcO - 
0.023 ± 0.004 cm"-'. As the maximum of EMpSin|fe| equals 
EMcSin|/7|, we have nliQ)h/2 = fonl^hll = 2.8 + 0.3cm-''pc 



^ Note that BMM defined h as the scale height of nl{z), thus their h 
equals our h^ = h/2. 



(see Sect. 2); hence nl(0) = 0.0060 + 0.001 lcm-^ WeO = 
0.26 + 0.03 cm-3 and /o = 0.09 + 0.02. 

We then analyzed the data in our sample using h^ - h/2 = 
0.47 kpc in eq. (|5]i, leaving out 4 pulsars in globular clusters at 
|z| > 3 kpc. As only a small part of their lines of sight passes 
through the electron layer, the mean electron densities along the 
line of sight are unrealistically low. Figure |2] shows the distri- 
bution of /iMp sin \b\(z) for all pulsars. The steady increase of 
EMp sin \b\ with |z| agrees well with the expected variation (full 
line), but the scatter is larger than in Fig.[T] 

4. Results 

In this section we investigate the dependencies of («e), («g), A^c 
and Fy on |z|, and the relationship between Fy and A^c- The sta- 
tistical treatment is the same as used by BMM and we refer to 
that paper for details. The results are given in Table|2] 
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y=21.7*(1-exp(-x/0.93)) 
y=21.1*(1-exp(-x/0.83)) 



|Z| (kpc) 

Fig. 1. Dependence of the dispersion measure perpendicular to 
the Galactic plane, DM sin \b\, on distance to the plane, \z\, of the 
final sample of 38 pulsars. Circles: pulsars with AD/D < 0.2; 
triangles: pulsars with 0.2 < AD/D < 0.5. Full line: the two- 
parameter fit to the data; dashed line: fit to the original sample 
of 78 pulsars, shown for comparison. 
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Fig. 2. Dependence of the corrected emission measure perpen- 
dicular to the Galactic plane, EMpSin\b\, on distance to the 
plane, |z|, of the final sample of 38 pulsars. Circles: pulsars with 
AD/D < 0.2; triangles: pulsai's with 0.2 < AD/D < 0.5. Full 
line: dependence expected for the scale height of 0.47 kpc used 
in eq. ©. 




DM (cm pc) 



Fig. 3. Comparison of corrected emission measure, EMp, and ob- 
served dispersion measure, DM. Circles: pulsars with AD/D < 
0.2; triangles: pulsars with 0.2 < AD/D < 0.5. Full line: bisector 
fit given in Table |2] 
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Fig. 4. Electron densities averaged along the line of sight and 
average volume filling factor as function of height above the 
Galactic plane, |z|. a) Average density («e) = DM sin \b\/\z\. The 
dashed line is the exponential fit to the data at \z\ < 1.2 kpc. 

b) Average of the square of the density (n^) = EMpSin|fe|/|z|. 

c) Mean density in clouds N^ - EMp /DM. d) Average volume 
filling factor Fy - DM^/(EMpD). Full lines indicate the expo- 
nential fits given in Table |2] Circles: pulsars with AD/D < 0.2; 
triangles: pulsars with 0.2 < AD/D < 0.5. 



We first show in Fig.|3]that EMp is well correlated with DM: 
the bisector fit yields EMp = (0.13 + 0.04)DM''-'**°°^ with very 
high significance. This indicates that both quantities probe the 
same ionized regions along the line of sight, although the beams 
used for the Ha observations (up to 1°) also sample regions 
located around the single sightline to the pulsar. Variations in 
electron density across the beam will contribute to the scatter in 
Fig. [3] The components of EMp and DM along |z| correlate less 
well than those along the Galactic plane (see Table |2]), because 
of the different scale heights of n^(z) and ne(z)- 



4.1. Dependence ofin^), (nl), Nc and Fy on height 

Figure |4] presents the various densities and Fy as a function 
of |z| in the logF - |z| plane. Although these variables repre- 
sent averages along the line of sight, we approximated their |z|- 
distributions by exponentials with scale height H. The fits are 
listed in Table |2] 

The distribution of {n^){z) in Fig. |4^ is effectively constant 
up to |z| ^ 1 kpc (dashed line), a well known fact first noted by 
Weisberg et al. ( |1980l l. Beyond this height (ng) decreases. A fit 
to all points (full line) gives a mean value at the midplane of 
("eX) - 0.019 + 0.02 cm"-', in fair agreement with the expected 
value («e){) = ne(0) = 0.023 +0.004 cm"^ derived in Sect. 3. The 
scale height is less well determined, though. The spread in the 



Elly M. Berkhuijsen and Peter Miiller: Densities and filling factors of the DIG in the Solar neighbourhood 
Table 2. Statistical relationships derived for the sample of 34 pulsars 



X 



Fit 



Corn 
coeff. 



DM 


EMp 


DMcos|fo| 


EMp cos \b\ 


DMsin|i| 


EMp sin \b\ 


D 


in.) 


\z\ 


(n.) 


\z\ 


inl) 


\z\ 


Nc 


kl 


Fy 


kl 


Lc sin \b\ 


A'c 


Fy 


A'c 


L, 


A'c 


Lc sin \b\ 



bis' 



lin^ 
exp' 



pow 



13+""^ 

^ -^-0.03 

13+""' 

" -^-0.02 

0.21 ±0.04 

0.018±0.002 

0.019±0.002 

0.0045±0.0007 

"•^^-0.03 

077+""2" 

"■"' ' -0.016 

80±20 

QQ,, +0.003 
"■"^^-0.002 

74+20 
2±1 



student 

t6 



1.15±0.07 


0.90 ± 0.08 


11.8 


1.15+0.05 


0.93 ± 0.06 


14.7 


0.98+0.09 


0.82 ±0.10 


8.3 


(2 ±6)10-'' 


0.52 ±0.15 


3.5 


-3- 


0.72 ±0.12 


5.9 


-i<i 


0.77 ±0.11 


6.9 


-2.6!i;? 


0.76 ±0.12 


6.6 


[+23]5 


0.63 ±0.14 


4.6 


1.06+0.15 


0.83 ±0.10 


8.5 


1.20±0.13 


0.88 ± 0.09 


10.2 


0.85±0.35 


0.75 ±0.12 


6.5 


-1.6±0.3 


0.82 ±0.10 


8.0 



Units: DM in cm"' pc, EM and EMp in cm"^ pc, {n^} and A'c in cm"', («^> in cm"*, 
kl in kpc, Le in pc 

1) Bisector fit to powerlaw Y = aX'' in log X - log Y plane. 

2) Linear fit; regression line of Y on X, where Y = a + bX. 

3) Exponential fit; regression line of In Y on kl with slope l/b, where b is the scale 
height H in kpc. 

4) Powerlaw fit; regression line of log Y on log A". 

5) Undetermined: errors > l/h. 

6) Student test: for a sample of 34 pulsars the correlation is significant at the 3cr level if t > 3.3. 
All fits are ordinary least-squares fits with errors of one standard deviation. 



data clearly decreases away from the plane, as was also noted by 
BMM. This is not due to longer pathlengths, because the spread 
in the distribution projected along the plane remains constant 
(not shown). It just shows that the variety in electron density is 
larger near the Galactic plane than away from the plane. 

The other panels in Fig. |4] show the dependencies of (nl), 
Nc and Fy on |z|. All quantities vary considerably about the fit- 
ted lines. The mean density in clouds slowly decreases from 
0.24 cm'^ at |z| = Okpc to 0.16cm~^ at \z\ - Ikpc, but at all 
heights values between 0.08 cm"^ and 0.5 cm"^ occur The mid- 
plane values of the exponential fits in Table |2] agree to within 
errors with those of nl{Q), n^o and /o derived in Sect. 3 and the 
scale height H of («g) is indeed about half that of («e) and A^c- 
This agreement shows that the data are internally consistent. 

The slight increase in Fy(z) in Fig.|4}l is statisically insignif- 
icant, thus Fy(z) is effectively a constant of about 0.08 + 0.02. 
This differs from the clear increase in Fy(z) found by BMM for 
a much larger sample of pulsars. The spread in Fy seems largest 
near the Galactic plane, but this impression is mainly due to three 
extreme points of nearby pulsars at low k|: two near Fy = 0.6 
and one near Fy = 0.07. At \z\ > 0.3 kpc the spread stays within 
a factor of 4. 

If the mean filling fraction remains about constant along the 
line of sight, the total pathlength through the ionized regions 
will increase linearly with distance. Lg sin \b\(z) indeed increases 
nearly linearly from about 7pc towards |z| =0.1 kpc to about 
80 pc towards |z| - 1 kpc (see Fig.|5]and Table|2]). 

4.2. Dependence of Fy on N^ 

One of the indications for turbulence in the ISM is an in- 
verse relationship between filling factor and mean cloud density. 
This is not only expected on theoretical grounds (Fleck [19961 
Elmegreen 1 1 9991 ), but has also been obtained from simulations 
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Fig. 5. Dependence of the total pathlength through ionized re- 
gions perpendicular to the Galactic plane, L^ sin \h\, on height 
above the plane. Circles: pulsars with ADjD < 0.2; triangles: 
pulsars with 0.2 < AD/D < 0.5. Full line: powerlaw fit given in 
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(Elmegreen 1 19971 Kowal & Lazarianl 2007l l. Observational evi- 
dence, however, is scarce (see Sect. 1). 

Since («e) is about constant (Fig.|4^) and (rig) = FyN^, we 
expect an inverse relationship between Fy and A^c. We present 
this relationship for our sample in Fig. |6] The correlation be- 
tween Fy and A^c is very good and indeed nearly inverse: the 
power-law fit yields Fy{Nc) = (0.011 + 0.003)A^,"'-2°^"'^ with a 
correlation coefficient of 0.88+0.09, and covers about 1.5 decade 
in N^ (0.05 - 1 cm"3) and Fy (0.4 - 0.01). We compare this re- 
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Fig. 6. Dependence of the average volume filling factor, Fy, on 
the mean density in clouds, A^c. for the sample of 34 pulsars in 
the solar neighbourhood. Circles: pulsars with AD/D < 0.2; tri- 
angles: pulsars with 0.2 < AD/D < 0.5. The full line shows the 
powerlaw fit given in Table |2] 



suit with earlier determinations and discuss its meaning for the 
density structure in Sect. 5.3. 

As Fy{Nc) is well defined, we looked for variations with lon- 
gitude, distance along the plane and height above the plane that 
may indicate variations in structure in the DIG. We only found 
a dependence on height: at \z\ < 0.3 kpc the relationship is con- 
siderably steeper and the filling factor for Nc - I cm"^ consid- 
erably smaller than at |z| > 0.3 kpc (see Fig. [T] and Table [3]). 
The two lines cover nearly the same range in Nc and cross near 
A^c = 0.15 cm"-', but at |z| > 0.3 kpc the spread in the data is 
larger than at lower \z\. On examination of the much larger sam- 
ple of BMM, it appears to show the same trend (see Table [3]). 
The large number of pulsars even permits a division of the lower 
|z|-interval yielding exponents of -1.6 + 0.2 for |z| < 0.2 kpc 
and -1.1 +0.1 for 0.2 < |z| < 0.3 kpc. It seems that in the 
thin Galactic disc Fy depends more strongly on A^c than above 
|z| = 0.2 kpc, with an exponent significantly smaller than -1. 
This indicates that the structure of the DIG in the thin disc dif- 
fers from that away from the plane. This may not be surprising. 
In the thin disc the activities of stellar winds, classical H II re- 
gions and supernova remnants largely determine the structure 
of the DIG. They may change or destroy the structure of the 
turbulence that is typical for more quiet regions away from the 
Galactic plane. Interestingly, Cordes & Lazio (2002) derived a 
higher fluctuation factor for the thin disc than for the thick disc 
in their model of the electron density distribution. There is also 
evidence from measurements of interstellar scintillation that dif- 
ferent types of turbulent spectra exist in different regions of the 
Galaxy (Shishov et al. i2003i) . 



5. Discussion 

5.1. Comparison with BMM 2006 

It is interesting to compare our results with those of BMM who 
took the distances to the pulsars from the NE2001 model of 
Cordes & Lazio (2002) and used an average correction for the 
absorption of the Ha emission that only depends on latitude. 
Applying their method to our sample, we recalculated the main 
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Fig. 7. Dependence of the relation Fy{Nc) on distance to the 
Galactic plane |z|. Left: |z| < 0.3 kpc; right: |z| > 0.3 kpc. The 
lines are powerlaw fits given in Table[3] 



relationships using a scale height of 0.47 kpc to correct for Ha 
emission from behind the pulsar. Table |4] shows that the result- 
ing relations are nearly identical to those derived in Sect. 4, but 
with larger errors. This means that the less accurate distances and 
absorption corrections used by BMM have not influenced their 
statistical results. Therefore, we may compare results from our 
small sample with those from the much larger sample analyzed 
by BMM with confidence. 

A striking difference between the two samples is the be- 
haviour of Fy(z). Where BMM found a significant increase of 
the filling factor to about 0.2 towards |z| = 1 kpc, Fy(z) derived 
from our sample remains constant within the errors at a value of 
about 0.08. This difference must be due to the different locations 
in the Galaxy of the pulsars in the two samples. For example, 17 
out of the 34 pulsars (50%) in our sample are within 1 kpc from 
the Sun, whereas only 44 out of 157 pulsars (28%) of the BMM 
sample are within this distance. Since the ISM is a highly vari- 
able medium, differences between samples are to be expected. 
Hydromagnetic simulations of the evolution of the ISM by de 
Avillez & Breitschwerdt ( I2005I I may illustrate this point. Their 
Fig. 2a shows the density variation in a cut of 1 kpc along and 
10 kpc perpendicular to the Galactic plane. Although the filling 
fraction of the DIG generally increases with height, there are 
also regions where it remains constant or even decreases away 
from the plane. 

5.2. Extent of ionized regions 

The total pathlength through ionized regions perpendicular to the 
Galactic plane increases nearly linearly from about 7 pc towards 
|z| - 100 pc to about 40 pc towards |z| = 500 pc and about 80 pc 
towards |z| = 1 kpc (Fig.|5]l. 

In the ISM dense regions are usually smaller than less dense 
regions, which also applies to the DIG (see Fig. [SJ. Although 
the distribution of Lg sin \b\ with A^c shows considerable spread, 
it is well fitted by a power law (see Table|2]i. Regions of density 
Nc = I cm"-' together occupy about 2 pc in the |z|-direction and 
tiiose of A^c = 0.05 cm"-^ about 300 pc. 

We calculated the same relationship for the BMM sample 
which gives a somewhat larger extent for the higher densities 
(see Fig. [8). This power law is close to the one indicated by 
the distribution of 192 pulsars presented by Hill et al. (I2007I I. 
However, these authors did not apply any correction to the ob- 
served emission measures, which leads to a flattening of the dis- 
tribution and a shift to higher values of A^c for the same path- 
length (see Fig. 5 in BMM). 
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Table 3. |z|-dependence of the relation Fy(Nc) - aN'^ 



Sample 



kl 

[kpc] 



N 



Corr. 
coeff. 



student 
t' 



This 


0.0- 


-0.3 


0.0070±0.0018 


-1.55±0.15 


15 


0.94±0.09 


10.2 


work 


0.3- 


-2.0 


0.021 ±0.006 


-0.80±0.17 


19 


0.80±0.14 


5.6 


BMM 


0.0- 


-0.2 


0.0078±0.0020 


-1.62+0.16 


11 


0.97±0.08 


12.7 




0.2- 


-0.3 


0.0183±0.0037 


-1.11±0.12 


24 


0.91±0.09 


10.1 




0.0- 


-0.3 


0.0137±0.0023 


-1.29±0.10 


35 


0.92±0.07 


13.9 




0.3- 


-2.0 


0.0187±0.0014 


-1.04±0.03 


122 


0.95±0.03 


31.9 



1) student test: for a sample oi N > 11 pulsars the correlation is significant at the 
3cr level if t > 4.0 



Table 4. Comparison with BMM for our sample of 34 pulsars 



Relation 


Our data 




BMM 






a 


b 


a 


h 


(nl) =fle-W'' 


0.0045±0.0007 


1 4+0.6 


0.0053±0.0011 


1 4+' 5 


Nciz) =ae-l~l/'' 


0.24±0.04 


2 6+-" 


0.26+0.06 


9+10 
^-1 


F,{z) =flel^i/'' 


0.077±0.018 


[23]' 


0.077+0.023 


[5]' 


F,{N,) = aN'^ 


0.011±0.003 


1.20±0.13 


0.014+0.003 


-1.18+0.11 



1) Undetermined: 1/b smaller than errors 
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Fig. 8. Dependence of the total pathlength through ionized re- 
gions perpendicular to the Galactic plane, Le sin \b\, on the mean 
density in clouds, Nc- Circles: pulsars with AD/D < 0.2; trian- 
gles: pulsars with 0.2 < AD/D < 0.5. Full line: powerlaw fit 
given in Table|2j dashed line: powerlaw fit to the sample of 157 
pulsars of BMM, L^ sin \b\iNc) = (5.5 + 0.8)A^-'-^^2±o.07^ 



5.3. The Fy-Nc relation and density structure 

In spite of our small sample of pulsars, the nearly inverse corre- 
lation between Fy and Nc shown in Fig.|6]is very tight, reaching 
a correlation coefficient of 0.88 + 0.09 (see Table|2]i. It does not 
change with distance along the plane and, as BMM explained, is 
insensitive to errors in emission measure. The near constancy of 
(we) = FyNc along the line of sight (see D - {ric) fit in Table|2]), 
first noted by Weisberg et al. ( 119801 ), was an early indication of 
an inverse relationship between Fy and A^c- 

We first compare our Fy-Nc relation with earlier determina- 
tions available in the literature, and then discuss what it tells us 
about the density structure in the DIG. 
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Fig. 9. Dependence of the volume filling factor, Fy, on mean 
density in clouds, A^c, as derived by different authors. Full line: 
this work; dashed line: dependence for the DIG within about 
3 kpc from the Sun (BMM); dotted line: dependence for the DIG 
in the bright emission ring in M31 (Berkhuijsen 2004); dash-dot 
line: dependence for clouds of diffuse dust within 400 pc from 
the Sun (Gaustad & Van Buren [1993]). The parameters of the 
lines are given in Table |5l 



5.3.1. Comparison with oilier data 

Figure|9|and Table|5]show our Fy-Nc relation together with ear- 
lier determinations. For the same mean density in clouds, we find 
somewhat smaller filling factors than BMM which agree with 
the lower points in their distribution. The smaller errors in the 
BMM relation are not only due to the larger sample, but also to 
the smooth distribution of (we) - DM/D of the electron density 
model of Cordes & Lazio (2002) used by BMM. Berkhuijsen 
(2004) obtained the Fy-Nc relation for the DIG in the galaxy 
M 3 1 from a comparison of rotation measures of the polarized 
continuum emission (giving (/ig) with magnetic field strengths 
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of Fletcher et al. 120041 ) and thermal radio emission (giving (n^)). 
These data refer to the DIG in the bright emission ring in M 3 1 . 
The good agreement between the M 3 1 and the MW results sug- 
gests that the inverse correlation between Fy and Nc is generally 
valid in the DIG in galaxies. 

Roshi & Anantharamaiah (120011) observed radio recombina- 
tion lines from the inner Galaxy at low latitudes. They derived 
filling factors of < 0.01 for extended regions of diffuse ionized 
gas with densities of 1 - 10 cm"^, which agree well with the ex- 
tension of the curves for the DIG in Fig.|9] 

Thus the inverse correlation between Fy and A^c at least holds 
for the density range 0.03 - 10 cm"^. Cordes et al. (1985) es- 
timated a filling factor of 10"'*''=^"-'' for clumps of about 1 pc 
size causing scattering of pulsar signals. Our data predict this 
filling factor for clumps of density Nc - 50 cm"^ and size 
Le/w = 0.9 pc, where m is the number of clumps on the line 
of sight. The relation of BMM gives clumps of Nc = 100 cm"^ 
for the same filling factor and sizes. So these clump properties 
also seem to follow the Fy-Nc relations derived for the DIG. 

Furthermore, Gaustad & Van Buren (T993S obtained a sim- 
ilar relationship for clouds of diffuse dust within 400 pc from 
the Sun, with mean densities of 1 - 10 cm"^. Figure |9] shows 
that these clouds have about 3 times higher filling factors than 
ionized clouds of the same mean density. Interestingly, their 
data agree very well with those of Pynzar (|1993l also shown in 
Berkhuii sen 1 1 998 J derived from recombination fines. 

5.3.2. Density structure of the DIG 

The general validity of the inverse Fy-Nc relation indicates that 
it describes a basic property of the DIG, possibly even of the en- 
tire ISM (Berkhuii sen 1 1 9991 ) . At least two mechanisms could be 
responsible: thermal pressure equilibrium and turbulence caus- 
ing a fractal density structure. 

Thermal pressure equilibrium leads to an inverse Fy-Nc re- 
lation if it were widespread, but it seems only locally valid. In 
the MW large fractions of gas are observed in unstable regimes 
and have also been found in simulations of a turbulent ISM. 
Thermal pressure equilibrium of clouds appears to be of mi- 
nor importance in the presence of turbulence (see Elmegreen & 
Scalo 2004 and references therein). 

Elmegreen ( 119981 [1999) discussed the properties of diffuse 
ionized gas in a pervasive fractal ISM. One of these properties 
is an inverse correlation between volume filling factor and gas 
density, which is expected on theoretical grounds (Fleck [19961 ) 
as well as from simulations (Elmegreen l 19971 Kowal & Lazarian 
120071 ). As this leads to a constant average density along the line 
of sight, the observed near constancy of (we) with Nc and Fy in- 
versely varying over 1.5 decade may indicate that the density 
distribution of the DIG is fractal. Please note that the near con- 
stancy of (rtc) alone should only be regarded as a first indication, 
because the fits in Table[2]show that the statistical significance of 
the near constancy of (ric) is much lower than that of the nearly 
inverse Fy-Nc relation. 

The fractal medium is characterized by filamentary, clumpy 
structures that tend to cluster together, with holes inside and 
large voids between them. The filling factor along the line of 
sight is determined by the outer regions of the filaments, i.e. by 
the largest scales, because the small, dense clumps in the inner 
parts hardly contribute. The voids occupy more space than the 
filaments, thus the filling factor of the filaments is small. The val- 
ues of Fv - 0.1 (Sect. 4.1) and Fy ^ 0.2 (BMM) obtained when 
looking perpendicular to the Galactic plane towards |z| = 1 kpc 
are consistent with this picture. 



The tendency of filaments to cluster into larger complexes 
reduces the number of fractal clouds, A^f , along the line of sight. 
Elmegreen (1998) estimated A^f - 3/kpc locally. For Fy - 0.1 
the total pathlength through the ionized regions is 100 pc/kpc, 
giving a line of sight through one fractal complex of about 30 pc 
and a mean density A^c of 0.15 cm"-' (from Fig. [6]). Using the 
results of BMM gives Fy = 0.2, a fractal complex of about 70 pc 
along the line of sight and Nc = O.lOcm"^. Interestingly, these 
pathlengths are in the range of 10-100 pc that Ohno & Shibata 
(1 1993 J estimated for the cells of ionized gas causing DM and 
rotation measure RM towards pairs of pulsars close on the sky. 

According to Elmegreen (119981 ), pervasive fractal structure 
only develops in regions outside the influence of SN shells, sites 
of star formation, chimneys etc. that are concentrated near the 
Galactic plane. This may explain our result of Sect. 4.2 that the 
exponent of the Fy-Nc relation is significantly smaller in the thin 
disc (|z| < 300 pc) than further from the plane where it is close 
to -1 (see Table[3]and Fig. [7]). 

We conclude that our results on the mean electron densities 
and volume filling factors of the DIG outside the thin disc are 
consistent with a fractal ionized medium caused by turbulence, 
while the structure in the thin disc is dominated by the effects 
of SN shocks, stellar winds, star formation and other forces. The 
difference in density structure between these two regimes could 
be further analyzed when more pulsars with measured distance 
become available. 

Recently, further observational evidence for turbulent struc- 
ture of the DIG and of the diffuse atomic gas has been presented. 
Hill et al. (120081 ) found that the probability distribution function 
(PDF) of the emission measures observed perpendicular to the 
Galactic plane at \b\ > 10° is lognormal as is expected for a tur- 
bulent medium from MHD simulations of the ISM (Elmegreen 
& Scalo 2004). Furthermore, Berkhuijsen & Fletcher (120081 ) 
showed that the PDFs of the quantities (ng), {n^}, Nc and Fy 
derived above as well as that of the mean H I density along the 
line of sight towards stars, {rim}, are lognormal, consistent with 
a turbulent origin of density structure in the diffuse gas. 

6. Summary 

We have used pulsar dispersion measures (Manchester et al. 
2005 ) and extinction-corrected emission measures (Finkbeiner 
2003; Dickinson et al. 2003) towards 38 pulsars with known 
distances for a statistical study of several parameters of the dif- 
fuse ionized gas (DIG) in the solar neighbourhood. The emis- 
sion measures were also corrected for contributions from be- 
yond the pulsar distance. To avoid regions with strong absorp- 
tion (> 1 magnitude) and contributions from HII regions most 
pulsars in the sample are at Galactic latitudes \b\ > 5°. The sta- 
tistical results are collected in Tables[2]and[3] 

Our main conclusions are: 

1 . From the scaling of dispersion measures with distance per- 
pendicular to the Galactic plane, we find a scale height of 
the ionized layer of 0.93 ±0.13 kpc and an electron density 
at the midplane of 0.023 + 0.004 cm"-', in good agreement 
with earlier determinations. 

2. Dispersion measure DM and corrected emission measure 
EMp are well correlated, indicating that they probe the same 
ionized regions. We may then combine them to derive the 
average densities along the line of sight, («e) and («g), the 
mean electron density in clouds in the line of sight, A^c, the 
volume filling factor of these clouds, Fy, and the total path- 
length through the ionized regions, Lg. 
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Table 5. Comparison of F^ - Nc relations Fy(Nc) - aN^ 



Sample 


a 


b 


Range A^c Range F, 
[cm-3] 


Ref. 


DIG 

MW N=34 
MW N=157 
M31 

DIFFUSE DUST 
D<400pc 


0.0114+0.0025 - 
0.0184+0.0011 - 
0.016+0.004 

0.056+0.020 


-1.20+0.13 
-1.07+0.03 
-1.24+0.30 

-1.25+0.40 


0.05-1 0.4 -0.011 
0.03-2 0.8 -0.009 
0.19-1 0.13-0.016 

1 - 10 0.06 - 0.003 


1 
2 

3 
4 


References. (l)Th 


is work ; (2) Berkhuijsen, Mitra, & 
04l (4) Gaustad & Van Buren|l993 


Mullerl2006l 




(3) Berkhuijsen|20 


1 





3. The total extent of the ionized regions perpendicular to the 
Galactic plane increases linearly reaching about 80 pc to- 
wards |z| - 1 kpc. 

4. The filling factor F^ remains essentially constant with |z| at 
a mean value of 0.08 + 0.02. Whether F^iz) is constant or 
systematically increases with |z| depends on the regions in 
the Galaxy probed by the pulsar sample. 

5. The average electron density («e) is about constant with a 
spread of a factor % 2 about the mean value of 0.018 + 
0.002 cm"^. Since {n^) - F^Nc, an inverse relationship be- 
tween Fy and A^c is expected. 

6. We derived the relation Fy{N,) = (0.01 1 + 0.003)A^"'-2°^°'^ 
which covers about 1.5 decade in A^c (0.05 - 1 cm" ) and F^ 
(0.4 - 0.01) (see Fig.|6]l. This is in good agreement with ear- 
lier determinations of the Fy-N^ relations for the DIG (see 
Fig.0. 

7. Near the Galactic plane the dependence of Fy on N^ is signif- 
icantly stronger than away from the plane. This indicates that 
the thin disc has a different turbulent structure than regions 
away from the plane. 

8. The inverse relationship between Fy and Nc, and hence also 
the near constancy of («e), are consistent with a fractal den- 
sity distribution in the DIG caused by turbulence which dom- 
inates the structure outside the thin disc. 
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